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NATIONAL ADVISORY COMMITTEE FOR AERONAUTICS. 

TECHNICAL MEMORANDUM NO. 46T8. . 

CHOICE OF PROFILE FOR THE WINGS OF AN AIRPLANE. 
By A. Toussaint and E. Carafoli. 

The choice of the profile for the wings of an airplane is 
a problem which should be solved by a scientific method based on 
data obtained by systematic experimentation. 

The problem, in its present form, may be stated as follows: 

"To find a profile which has certain required aerodynamic 
characteristics and which encloses the spars, whose number, 
dimensions and separating distance are likewise determined by 
structural considerations." 

These conditions, imposed in the choice of the profile, 
result from the specifications to be satisfied in the static 
tests. At present, the static test, corresponding to the case 
of accelerated flight at limited speed, requires the knowledge 
of the moment of the aerodynamic resultant at the angle of zero 
lift, and the possibility of controlling the magnitude of the 
corresponding absolute coefficient within more or less extens- 
ive limits. Let us add also that, according to practical con- 
structional procedure, . this absolute coefficient (0m o ) should 
be very small. 



*"Le choix du profile des ailes sustentatrices , " from L 1 Aero- 
naut ique, December, 1927, 
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Moreover, .a knowledge of the pressure distribution around 
the profile and along the wing span is necessary in all cases' 
of static tests, since the diagrams of the load distribution of 
the static tests must be identical with the diagrams of the, 
aerodynamic pressures. 

Lastly, the profile should also have a profile drag as 
small as possible and good lifting qualities. 

The search for a profile answering the various requirements 
can sometimes be made by consulting the collections of profiles 
investigated in the aerodynamic laboratories. At present, how- 
ever, these data are generally insufficient, because they do not 
give the aerodynamic-pressure distribution and also because the 
investigated profiles have too large an absolute coefficient 
°m 0 > or dimensions and forms which do not satisfy the struc- 
tural requirements. 

In this case, the constructing engineer had, in the past, 
no other resource than that of adopting one or more profiles 
which he drew empirically and which he subsequently caused to be 
investigated in the aerodynamic laboratory. This long and bur- 
densome process did not always yield the desired solution. Be- 
sides, the process had to be repeated very frequently without 
ever leading to general and conclusive results. 

The only rational solution of the search for an airplane 
wing profile is found by scientific methods of drawing and cal- 
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collating so-called "theoretical" profiles. It is thus possible 
to find an infinite number of profile forms which answer all the 
proposed problems. For all these forms it is possible by start- 
ing from the data of the diagram to calculate the following es- 
sential characteristics: 

1. The variation of the lift coefficient C 2 , as a func- 
tion of the angle of attack, and, particularly, the value of the 
angle of zero lift. 

2. The variation of the moment coefficient C m jv, of the 
aerodynamic resultant with respect to the leading edge, and, 
particularly, the value C mQ of this moment for the angle of 
zero lift. 

3. The aerodynamic pressure distribution along the wing 
chord for all angles of attack used in practice. 

TLstneory for wings or cells of limited span allows, the ex- 
tension of these characteristics to all sections of the wing 
along the span. 

Conversely, this method renders it possible to find the 
profiles answering the required conditions, which is generally 
the" Sase in practice, as regards the value of the moment coeffi- 
cient 0m o # 

Convinced of the practical utility of these methods of draw- 
ing and calculating theoretical profiles, we have written quite 
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an exhaustive treatise on the subject.* At the request of the 
editor of this magazine, we have consented to give, with the ex- 
ception of too greatly elaborated theoretical considerations, 
the essential elements for drawing theoretical profiles to an- 
swer actual problems. The same opportune question has been fre- 
quently put to us by engineers. We have had the satisfaction of 
learning from those interested, that the wing profiles, thus 
conceived and investigated, have verified the theoretical pre- 
dictions in laboratory and flight tests. 

Elementary Principles for Drawing a Theoretical Wing Profile 

A theoretical profile is derived by conformal transformation 
from a circular contour. 

In order to effect this conformal transformation of the cir- 
cle into a profile a transformation function is used, which must 
satisfy certain mathematical conditions. The corresponding cal- 
culations are made by the use of complex variables. 

Let the plane of the generating circle with center M be 
Oil) (rigc, l). Each point P 1 , with coordinates % and T] 
on this plane, will be represented by the complex variable:^ 

C = t + ir\ = p e 1 ^, with P = 0P» (l) 

*A. Toussaint and E. Carafoli, ,.«Theorie et traces des profils 
d* ailes sustentatrices, " in course of publication in the "Bulle- 
tin de la Ohambre Syndicale des Industries aeronautiques, Vol. V, 
Bulletins Nos. 1-3, have already appeared* 
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Let us now consider the plane 0 xy. Each point P, with 

coordinates x and y in this plane, will be represented by the 
complex variable: 

z = x + iy = r e 1<p , with r = OP- (3). 

The most general transformations function for accomplishing 
the conformal transformation of the circle in the plane £ into 
a profile in the plane z can be expressed as 

z = t + J 1 " + + •••• + • (3) 

The form and aerodynamic characteristics of the transformed 
profile depend on the number of terms of this function and on 
the parameters oc 1 , x 8 , x 3 .... x n which, in the general case, 
are complex quantities. In particular, the first parameter x x 
may be written 

oc x = c s e 1 s7 (4) 

and the half-amplitude (demi-argument ) 7, characterizes the 
direction, of the second axis of the profile, that is, the direc- 
tion of the relative wind for which the aerodynamic resultant 
passes through the center M of the generating circle, when the 
£ and z planes are superposed and the axes .0 £rj and 0 xy 
coincide. 

In practice, the parameters x 1 , x 3 , .. x n are determined 
in the following manner. 
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By assuming the roots (-X), X x , \ 2 . . . X n which nullify 
the derivatives || of the transforming function, $e obtain 
the equation 

= (i + f)(i-^)(i-^).... (1 _M) 

from which, by identification of the two developments, we obtain 

n 

X x + \ 3 + ... + >fc = S Xi = X (6) 

jc x = c* e 2i<Y . X» - £ Xi Xj (7) " 



n n 
2x 3 = — X Z Xi X-j + Z Xi X-j X]j 



(8) 



The roots (-X) X x , X s , , X n , in the general case, 

are complex quantities. It is customary, however, to take the 
root (-X) along the axis 0 I of the real quantities, its 
magnitude being arbitrarily chosen. This magnitude character- 
izes the profile dimension, that is, the scale of the drawing. 

Moreover, the roots must satisfy only the following theoret- 
ical conditions: 

1. The corresponding representative points L, L x , L 2 , 
L n must be within the circle. However, in the case of an air- 
foil terminating in a sharp trailing edge (which is the general 
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case in actual practice), one of the roots and only one is lo- 
cated on the circle, (incidentally, it may be remarked that it 
is possible to draw a profile without a sharp trailing edge or 
with only a more or less rounded edge. ) It is customary to choose 
the real root £ = -X, so that the corresponding representative 
point L is coincident with the point B» of the generating 
circle (Fig. l). 

2. The geometric resultant of the vectors 0L ly 0L 2 , .... 
0L n corresponding to the roots X 1 , \ a , .... \ n must be equal 
to OL' = X by virtue of equation (6). 

These theoretical restrictions of the choice of roots are 
supplemented by lessons taught us by the practice of drawing. 
Thus all roots such as L ± (Fig. 2) near the circle at A 1 will 
produce a slightly rounded point at" A on the profile. Simi- 
larly, a root such as L-j will produce a boss at J on the 
prof i le a 

Practical Method for Drawing a Theoretical Wing Profile 
with a Sharp Trailing Edge 

The practical method for drawing a theoretical Tiring profile 
with a sharp trailing edge, according to the preceding general 
principles, comprises the following operations: 

1. Draw the coordinate axes for the planes £ and z in 
coincidence, and let 0 xi and 0 yr) be these axes (Fig. .3). 
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2. Choose the characteristic magnitude X according to the 
dimensions of the profile to he obtained, knowing that the chord 
I of the profile will be approximately 4 X . 

3. Lay off OB' = X on the 0-x, %■ axis in the negative 
direction. The point B' is located on the generating circle 
and represents the real root £ = - \ . 

4. Choose the other roots X x , X 2 , X n (equal in 
number to the terms desired for the transformation function) so 
as to satisfy the abovementioned conditions. 

The center M of the generating circle is not yet known, 
but it is, in practice, not far from the origin 0, so that the 
condition relative to the position in the circle for the repre- 
sentative points L x , L 2 , L n can be realized with an ap- 
proximation, which in general is proved sufficient by the con- . 
tinuation. of the drawing. The parameters x lt x 2 , x 3 , x^ 
are then calculated by means of equation (5). 

5. Draw the second axis of the profile. It is known that 
the amplitude of the second axis is equal to half the amplitude 
of the parameter: 

n 

x x = c2 e l3<Y = X s - Z X± Xj (7) 
To that effect, we determine the point C representative 

x + \l ^ (7.) 
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since 0B» = - X, it will suffice to draw through the point 

B» a straight line B«C making with 0 x an angle equal to 

n 

the amplitude of the complex quantity Z X ± Xj . On this axis 
and with its sign we lay off a length B»0« such that 



n 

mod Z Xi Xj 



n 



Z Xi Xj 



A, A 

We then have 

c 2 = x± = X x 0C» , (9) 

and the second axis is the bisector of the angle 

B« OC =87. 

The quantity cs = x x , which appears in the expression 
for the moment of the resultant, can then be calculated by meas- 
uring 00' directly on the diagram, or by solving the triangle 
B» 00«. 

6. Draw the first profile axis for the purpose of obtaining 
a fixed moment coefficient 0 mo . 

The direction of the first axis corresponds to that of zero 

lift. 

We shall see farther on,, in connection with the aerodynamic 
characteristics of the profile, that the value of C mo is given 
by the expression 

Cfco * 8 it S£ (p -.?) (10) 
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in which I is the profile chord, the value of which is approx- 
imately 4 X, or better still 2 (X 4 00» ), and p is the ampli- 
tude of the first axis or axis of zero lift. 

For the required value of 0 mQ , we therefore have 

Following laboratory practice, we consider C mQ positive, 
when it tends to diminish the angle of attack. Thus, for 
Cm 0 > 0, we have p > 7; for C mo = 0, we have p = 7; and 
lastly, for C mQ < 0, we have P < 7. 

The first profile axis is the straight line B» MI drawn 
through the point B 1 and making with the axis Ox, I the angle 
P determined by formula (11). 

7. Draw the generating circle. The center M of the gener- 
ating circle is found on the first axis, which meets the axis 
Oy, f] at the point M 0 . The distance M 0 M is arbitrarily cho- 
sen very small. (This distance characterizes, in part, the rela- 
tive maximum profile thickness.) One may thus define the center 
M by adopting an arbitrary value, a little larger than unity, 
for the ratio 

v _ B«M a - • ■ 
* " ~T = X 



a = B'M is the radius of the generating circle. 

Lastly, we may characterize the position of the center 
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by the amplitude 6 of the direction OM. From the diagram 
(or by calculation of triangle B'OM), the value of the modulus 
= OM. 

The generating circle is drawn through the point B 1 on 
the axis 0-x, % and we may verify the required condition for 
the interior position of the points L x , L 2 , L n . 

8. Draw the auxiliary circle corresponding to the term 

x 

-± of the transformation function. 

In drawing the profile and the generating circle on the 
axes Ox' and Oy* so that Ox' coincides with the second axis, 
the position of any point P' on the circle and of the corre- 
sponding point P on the profile are characterized by the com- 
plex variables £' and z' respectively, so that we have 

z = z« e i7 and £ = V e i7 



from which, we derive 



■ i _ 



V + a -, a + .... (13) 



x 2 e~ aiy . x a e-* i7 



In this form, it is obvious that, in order to find a point 
P of the profile, it is first necessary to evaluate the sum 

It has been demonstrated that the auxiliary transformation 
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leads to another ■ circle (called auxiliary circle) defined as fol- 
lows. The center M x of this circle is located on the line 
0M X symmetrical to OM with respect to the new axis oy r at a 
distance 0M X calculable "by the formula 

Lastly, the radius a x of this circle is given by the 
equation 




(14) 



It can be shown finally that the auxiliary circle passes 
through the point 0' , so that 

a x = M^O 1 " 

9. In order to accomplish the partial transformation corre- 
sponding to the first two terms of the transformation function, 
it is sufficient to apply to the circles M and M x the con- 
struction of Trefftz, commonly applied in the case of Joukowski 
profiles, or profiles with two terms* 

Fox the sake of completeness, we shall briefly review said 
construction. 

With the origin 0 as center, (Fig. 4) a circle is drawn 
which is divided into equal arcs (generally of 30° or 15°) 
starting from the second axis Ox'. The radius vectors, corre- 
sponding to these arcs, define, on the circle M, the points 
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P'o, P'j , P'k, and, on the circle 1£ 1} the points Q'q,^. 

Q'lo which are conjugated two "by two for amplitudes of equal 
magnitude but of opposite signs. 

The parallelogram constructed upon two conjugate radius 
vectors, 0P' k and 0Q' k defines a point P k of the desired 
transformation. The construction on the axes Ox 1 and Oy 1 
leads to the addition or subtraction of the conjugate radius 
vectors. 

For example, 0P o = 0P' o + OQ'o and °P 3 = 0P' 3 - 0Q» 3 . 
We thus obtain just as many points P k as we wish, according to 
the magnitude of the adopted angular equidistance. 

10, To make the partial transformations corresponding to 
the other terms of the transformation function. The other terms 
of the transformation function have the general form 

according to the expression 12 in the system of coordinate 
axes Ox 1 and Oy' . According to equation (5), the parameter 
x n can be calculated in terms of the roots of |~| and may be 
written 

*n = on e 10 * 

Every point P' k of the generating circle is defined like- 
wise, with respect to the axes Ox 1 and Oy 1 , "by the expression 

£' = p k e i9 k 
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from which it follows that the amplitude <p n of the term is 

' ■— •- qp n = a n - (n + 1)7 - n 6 k (15) 



and the modulus is 



x n e " 



>i(n+x )y 



n 



cn 

' n 



x 



lfnl_ 
„ n 



(16) 



We can calculate <p n and -S for the various points 

P' 0 , P^ , P' k , previously determined according to the 

adopted angular equidistance. 

The partial transformation corresponding to one of the other 
terms of the function will then consist of a positional correc- 
tion, the amplitude of which, with respect to the axis Ox' is 
q^, or (cp n + 7) with respect to the axis Ox, the magnitude 
and sign of which are given by 



x 



n! 



>n 



Thus, for all points P k , corresponding to the first two 
terms of the transformation, the third term 



x 2 e 



~si7 

2 



V 

defines a positional correction P k P k (Pig. 5), the amplitude 
of which is I 
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wjth respect to Ox 1 , or 

<P 2 + 7 = cr - 27 - 26 k 

with respect to Ox. 

a, being the amplitude of x 3 calculable by equation (8), 
and 9 k being equal to k times the adopted angular equidis- 
tance, since it concerns a point P k , the antecedent of which 
is P« k on the generating circle, the magnitude or modulus of 
the correction is 

■x 5 ^ x a 

e? = (op k ) s * 

\x I 

Remarks .- In practice, we compute the value of ' ^ for 

P k 

a few points P' 0 , P^ , p« k and draw a curve of the varia- 

tions of this quantity, which is generally small. The other val- 
ues are then derived by interpolation, according to the graph- 
obtained. 

We can therefore easily make the various positional correc- 
tions corresponding to the partial transformations due to the 
terms of the function other than the first two. 

Remarks on the Practical Choice of the Number of Terms 
of the Transformation Function. 

In practice, it is not necessary to resort to a very com- 
plicated transformation function. One of three terms, such as 

x x*^ 

z = t + -f- ■ + -73 ( 32 ) 
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is capable of .furnishing an infinite number of profile forms 
with the required 0 aiQ , In this case, the drawing will reduce 
itself to the graphic construction for the first two terms 

v £' J 

and to a single positional correction due to the term 

— *-% 

However, experience shows that, in this case, the position 
corrections are very important. It follows that the envelope 
of the points corresponding to the first two terms ( V + 
is very different from the form of the final profile. It is im- 
portant, from the practical point of view, after the first con- 
struction ( £' + 2Jty to be sure of the form and approximate 
proportions of the final profile. We thus find, without draw- 
ings or useless computations, that the, profile sought will ef- 
fectively and very closely envelop the spars necessary for the 
construction. 

Under these conditions, we are frequently led to use a 
transformation function of three terms corresponding to a par- 
ticular case of a four-term function. This particular case 
consists in nullifying the third term in x s . The function is 
then written 

2 = t + y + o + |i (33) 
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In order to obtain this particular four-term function, it 
suffices to choose the roots \ lt X 2 and X 3 in such a manner 
that jc 2 = 0. 

We have already found that 

n n 

2x 3 = - X S \ ± X-j + 2 \ ± X-j \ k (8) 
When n = 3, the annulment of x 2 leads to the condition 

3 3 

X I \i \j = E \i \j X^, 



that is to say, 

MX X X 2 + X x X 3 + \ 3 \ 3 ) = \ a \ 2 x a . 

This condition is verified for 

\ = \ with X 2 = - X 3 (34) 

The roots are therefore so chosen that the corresponding 
representative points satisfy the conditions 

0L X = X directed along Ox 
and 0L 3 = - 0L 3 in any direction (Fig. 6). 

By this particular choice of roots, drawing the profile is 
rendered very simple by a graphical construction ( £' + <j*£\ and 
by a positional correction, generally very small and correspond- 
ing to the term 

x 3 a- 1 *** 



I 3 
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Remarks .- In practice, we can do without defining the roots 
.^2 = - If we assume an arbitrary value for the segment 

B'O 1 (Fig. 7), the parameters x 1 and x 3 of function 33 are 
determined. In fact, we have 

x x = \ (OC ) 

On the other hand, 

— x : x~ • 

Since 3x 3 = \ \ x \ a \ a = - \ 2 \ a a, 

we obtain _ y 3 

x 3 = + £ \a X(B'C ) = 3- (B'C ). 

The above exposition corresponds to the method of Von Mises. 
Later we shall give a simpler and more general construction, 
which will enable us to predict the evolution of the profile 
form. 

Example of a Profile Drawing 
with 0 mo Fixed in Advance 

We shall give (Fig. 8) the example of drawing a theoretical 
profile for which the condition 

0m o = 0.055 

is. required. We shall select for this drawing the three-term 
■ functions 

z =£+£i. + 0+ |§. (33) 
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Following the order previously indicated, we shall have to 

perform the following operations: 

1st. Drawing the coordinate axes Ox and Oy or Oi and 0*1. 

2d. Choosing \ We shall take X = 6 cm, in order to 

obtain a profile chord about 2* cm 
long. 

3d. We lay out OB' = - \ = - 6 cm on the negative di- 
rection of Ox. 

4th. Choosing the roots and computing x x and x 3 . 

According to what has preceded, we shall have 

\ 1 = \ = 6 cm 

and we assume that 

\ a = - \ 3 = 0.345 \ e 155 °= 2.075 e 15S ° 

The representative points L s and L 3 will be on the line 
of amplitude 55° at the distances 

0L 2 = - 0L 3 = 20.75 cm. 

From this choice of the roots, we obtain the following 
equation. 

ft = i - - = (i + -|) (i - (i - (i - Y) 



i Jl 0.119 X 2 e 111Q , 0.119 X? 



4- e ino° 



from which we deduce 
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jc x = X s (1 + 0.119 e ill0 °) = 36 (1 + 0.119 e ill0 '° ) 



x 3 



0.119 y e 1 ^ = _ 51 . 6 e i^oo 



5th. Drawing the second axis. We construct 

- *i = . 6 - 0.716 e illo ° 
In the direction of amplitude +110° we lay out 



B'C = - 0.716 cm, thus finding 

B'OC 6.6° 



= 3.3°, 



from which we get the second axis, passing through the origin 5. 

6th. Drawing the first axis. The amplitude B of this 
first axis will "be determined by the condition relative to 



0m o = 0*055. We lay out on the drawing, 00' = 5.77 cm, and 
calculate 

a y _ n n«w x (6 + 5.17) 2 x 57.3 _ p o 

° °' 055 * 6.28 x 6 X 5.77 " 2 ' 

We therefore have B = 7 + 2° = 5.3°, whence the first axis 
of amplitude 5.3° passing through B 1 . 

7th. Drawing the generating circle. The center M 1 is 
chosen in the direction OM with amplitude 6 = 52° + 7 , and 
we lay out on the drawing 

OM = |p | = 0.73 cm, 
bTM = a = 6.45 « 
0 2 = X 0G« = 34.6 cm. 
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8th. Drawing the auxiliary circle. Draw the line 0M Z 
symmetrically with OM with respect to Cy* (perpendicular to . 
the second axis). Then compute 

OMi = = ? 4,6 x 0.73 = 0.618 cm 

1 11 (6.45 f - 0.73 

from which we obtain the circle with the center Mj. and passing 
through C. 

9th. Partial transformation corresponding to the first two 
terms. Having adopted the angular equidistance of 15° we pro- 
ceed with the graphic construction! on the conjugate points of 
the circles M and M r . 

10th. Positional corrections corresponding to the term 
-£§. By computing the amplitude and the correction modulus for 
each of the points 1« , 3» , 3» , ... of the generating circle, 
we obtain 

(cp 3 +7) = no 0 - 3X3.3° - 3xn8 
with respect to the direction Ox, and 



• 1*3 1 . 51.6 . 
?S Pn 

in which 110° is the amplitude of x 3 and n6 = n x 15° is 
the amplitude of the point n 1 on the generating circle. 

Through each point resulting from the graphic construction 
on the circles M and M x the direction of amplitude (<P 3 +7) 



x 

is drawn, and the length — . is laid out in the negative direc 
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t< 

tion, since the corrections have the minus sign. As easily 
seen., the positional corrections are relatively small, so that 
the final profile differs very little from the envelope of the , 
points resulting from the first two terms. 

Figures 9 and 10 give two other examples of drawing theo- 
retical profiles. Figure 6 corresponds to the case of a three- 
term function of the form 33, while Figure 7 corresponds to the 
case of a three-term function of the form 32. In the latter case, 
the positional corrections are important and the envelope of the 
points, corresponding to the partial transformation due to the 
first two terms, furnishes no information whatever regarding 
the form of the final profile. 

To be followed by Technical Memorandum Ho. 469, containing 
the translation of the remainder of this article. 

Translation by 

National Advisory Committee 

for Aeronautics. 
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